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ON THE STRUCTURE OP THE UPPER ATMOSPHERE OP MARS 
ACCORDING TO DATA PROM EXPERIMENTS ON THE 
"VIKING" SPACE VEHICLES 


M. N. Izakov 


Altitude profiles of the concentrations 
of the atmospheric components measured by 
the on-board mass spectrometers during the 
descent of the "Viking" lander can be satis- 
factorily described by assuming that tempera- 
ture has a smoother profile, and the eddy 
mixing coefficients are smaller at altitudes 
of 120 to 170 km than those determined by 
the authors of the experiments in [1, 2] . 

The eddy mixing coefficient in the Mar- 
tian thermosphere during the "Viking" experi- 
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ments was equal to 5*10 cm /sec according to 
our interpretation. 


1. During the descent of the; lander of the "Viking-1-2" 
space vehicles into the Martian atmosphere, excellent measure- 
ments were made by means of mass spectrometers. These data 
were published by Nier and McElroy and their colleagues 
[1, 2] together with the data from measurements in the lower 
atmosphere [3, 4], This information constitutes a major con- 
tribution to our knowledge of the structure of the Martian at- 
mosphere. Using these data, investigations 'of the composition 
of the atmosphere of Mars and the processes occurring in it 
have been undertaken [5, .6] . 


However, in our opinion, certain details in the process- 
ing of the data discussed in [2] require further analysis, 
and certain conclusions are controversial. It is in connec- 
tion with these matters that the present review was undertaken 
E5, 6]. 

*Numbers in margin indicate pagination of original foreign text. 
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2. The altitude profiles of the temperature T(h) , 
which were computed In [2] from the altitude profile of the 
concentration of carbon dioxide gas 'fl > measured 

by a mass spectrometer, have a meandering form which, in the 
opinion of the authors, reflects the actual variations of the 
temperature associated with the passage of tidal waves in the 
atmosphere. To some extend, this can be the case; however, 
this can also be due in part to errors in measuring the con- 
centrations and calculating the temperature. 

The error in the temperature determined from the concen- 
tration profile is 60° at an altitude of 200 km, accord 

ing to the authors' estimate [2]. Then it quickly decreases 
with the altitude, and at altitudes from 140 to 160 km, Oi 
+ 5 - 10°, and then it again increases to 20 - 30° at 120 

km altitude. To achieve such a high degree of accuracy in the 
measurement of temperature at altitudes of 140 - 160 km, it is 
necessary to measure the concentrations with a precision of 
u’Ags 3-5 %. Is this possible in flight? In [2] it is 
assumed that it is possible to determine from various 

concentrations of a single component, determined from various 
mas:s peaks (for COg — 44, 22, 12). Here it is implicitly 
assumed that, as compared to calibration error, it Is possible 
to neglect: errors due to instability in' the characteristics 

of the apparatus in flight (in particular, the source of ions); 
errors due to calculating the concentrations In the atmosphere 
from the concentration in the equipment (the uncontrollable 
variations in the angle of attack may contribute to this error) 
and for components falling on some peaks, such as and CO, 

errors due to supplementary calculations to isolate them. But 
if all these errors are actually negligible, then, such a strong 
non-monotonic dependence of errors on altitude is not under- 
standable', A supplementary analysis of these' questions, which 
only the authors of the experiment can carry out, would be 
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extraordinarily valuable. It Is highly Improbable that the 
error In measuring the concentrations In flight can be markedly 
less than 25 - 30 % . 

When calculating the temperature from the concentration, 
using the barometric formula, was done In [2] , the temperature 
is determined essentially as the logarithmic derivative of the 
concentration, but during numerical differentiation of an 
empirical curve containing an error, the error in the result 
may be increased. Prom this point of view, it is more advis- 
able to integrate the concentration curve numerically, obtain^ 
ing the pressure, as was done in B] , and to obtain the temper- 
ature from the pressure. When a sufficiently precise computa- 
tion is performed by this method, the relative error in the 
temperature may remain the sa.me as the relative error in the - 
concentration. 

When differentiating as well as when integrating the pro- 
file of ji ./['ll , it is necessary to interpolate between the 
experimental points, since they are far apart. Knowing the 
error in determining the concentration 0 , it is advisable 

to construct a smoothed approximation curve (e.g., by the method 
of least squares), which deviates from the experimental points 
by not more than the magnitude of the error, rather than to use 
interpolation. This can help to exclude artificial waves from 

t 

the curve 

The profile of ^COp^ » com P uted with constant tem- 

perature relative to the height, deviates from the experimental 
points by not more than 30 % (Pig. 1), except for the lower joints 
at 130 km, which may have been influenced by the "satura- 

tion pressure" in the apparatus'”^] , and the upper point at 199 
km relative to "Viking-2." Therefore, we shall include the fact 
that at the time of the "Viking-l" experiment, the temperature 


of the Martian thermosphere was 180^40 ■ a > and at the 

time of the "Viking-2" experiment — 150^50 K. 

«y 

If the structure of the profiles of Kl^( ja ), ; T£ (t.) 
actually reflects the Influence of the atmospheric waves, 
then the smoothing curve should be close to the equilibrium 
profile in the absence of waves — the most important charac- 
teristic of the atmosphere, which must be known first of all. 

3. In [2] the very large values, K = (1 - 4). 10^ cm 2 /sec 
at an altitude of 170 km (Pig. 2) are obtained from the pro- 
files of the argon and nitrogen concentrations, > 

Yl ( 1/1 ) > and the temperature K) calculated from the 

profile of H. co £h ) 

** % 

In our opinion, the values of K contain a very large 
error, and at altitudes of 140 - 170 km, this error drastic- 
ally increases for the following reasons. In [2] K is found 
by means of the following formula: 

afriri* , 'a&nT 
•ah. + 

. which is the barometric formula in differential form for the 
transition layer between the homosphere and the heterosphere. 
Here 

& K! ( __ 1 

where is the altitude scale for the pressure of the 

olth component ; is the altitude scale for the 

pressure of bhe ^ th 'component during diffusion equilibrium 

in the heterosphere; jj _ is the mean altitude scale 
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during diffusion equilibrium in the homosphere . Clearly, 

***** 

in the homosphore where > and in the heterosphere 

where ^ , formula (1) is converted into the usual 

barometric formula for H or , respectively. Formula (1) 

is obtained from the continuity equation for a long-lived 
atmospheric component (i.e., a component for which the charac- 
teristic time of chemical reactions is much greater than the 
characteristic time of diffusion), if it is assumed that it 
is in a state of stationary distribution and if we neglect 
thermal diffusion (cf., e.g., [73). However, if acoustical- 
gravitational waves are propagated in the atmosphere, then 

the distribution deviates from the equilibrium distributions, 

Yl > 

and It is impossible to neglect the time derivative 

. *0 £ 
in the continuity equation, and in the case of large tempera- 
ture gradients, it is probably Impossible to neglect the 
thermal diffusion term also, and therefore formula (1) is not 
applicable. Nevertheless, let us assume that the measured 
profiles of T( n.) are close to the equilibrium dis- 

tributions. Then (and only then) is formula (1) valid. How- 
ever, here K as found from (1) has a large error.- Indeed, if 

we solve equation (1) for the desired K, we obtain 

i i 


K = 3, 


I jo 
HvL 


U 

rk 


X- - J- 
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Various authors differ in their use of the term "altitude 
scale." It is reasonable to define it as the altitude interval 
over which some quantity undergoes an Q -fold change, i.e., 

1 d&l 

rr = 7rr~Tu , • Here it is obvious that the altitude scales 
Hup .lY a*!' • .j. ,r 

for pressure and concentration are different: -j-r 

rls> Hi-, ’»■' * *• 

However, often, when speaking of an altitude scale, we have in 
mind an altitude sjjale for pressuring during diffusion equili- ’ 

brium -« 
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We see that K Is defined by a fraction whose numerator and 
denominator contain differences, each of which is encumbered 
with a significant error (since they contain the logarithmic 
derivative of the profile > the temperature cal- 

culated from the profile, and* the logarithmic derivative of 
the temperature); moreover, the altitude scales occurring in 
(2) do not drastically differ in value since the state is 
approximately the equilibrium state, and the region under 
consideration is not far from the homopause (where, by defini- 
tion, K = '“£j c » ) . It is completely obvious that the resulting 

error in the determination of K may be very large. But that 
is not all: to the error due to the structure of the expres- 

sion for K and the errors contained in its parameters is added 
the error resulting from the imprecise value of the mean 


molecular weight 


I';-' 
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which is contained in H. In- 


deed, in the range of altitudes under consideration, atomic 
oxygen becomes one of the basic components of the Martian at- 
mosphere whose concentration is not reliably measured by the 
mass spectrometer, but is- determined (with a certain error) 
from models. According to the models [6, 8], the concentration 

f 1 q is compared with the concentrations at 180 - 200 km 

CO- 


altitudes . 


We may assume that the errors mentioned mean that in the 
calculations presented in [2], the values of "drag after 

them" the values of K: indeed, it is clear from Pig. 2 that 

the function K(h) at altitudes greater than 130 km is very 
similar to the well-known function 

There is one further important circumstance. It is evident 
from (1) that if the molecular weight of some component M* 
is equal to the average molecular weight M, and therefore 
- H, then the equation degenerates, and it is impossible to de- 
termine K from it. But if is close to M, then again it 


is practically impossible to determine K, since very large 
differences in K yield very small differences in Pt 
Such is precisely the case for argon, since p'j^ 


CA,' 

= 40 is 


r**' 


close to M » 43 . 5 . Therefore, it is advisable to find K 
from the curve : 1^ is better to use this' curve 

to find T( k) 1 it is true that the error in the homosphere 
is somewhat larger than with respect to due to 

the greater difference between the molecular weight and the 
mean molecular weight. 


Summarising, we can say that in order to obtain a reli- 
able determination of the eddy mixing coefficient K by the 
method under consideration it is necessary, first of all, to 
use the profile of the concentration of a component with 
molecular weight which differs markedly from the mean (the 
profile of helium would be ideal) and, secondly, to have simul- 
taneously reliable measurements of atomic oxygen. 


4. To test the ideas discussed above we made such a cal- 
culation. In accordance with [1], we took the values of the 
concentrations of the components ilj at some initial altitude 
pi near the lower boundary of the region of the measurements 
and Pl,(h) were computed from the barometric formula, ex- 
pressed in the following form: 



C h 


Here the coefficients of diffusion M(-) GOo) were calculated 

in accordance with the data of [9], where considerable material 
from laboratory measurements is generalized, and. the function 

Aip -Aatp I ' ^ ft (wher ” constants for 

the given pair of molecules, cf. the table) is proposed. 





Parameters for < 
of diffusion of 
formula : 


CO - co 2 

0 2 — vj 0 2 

i\ r 2 - C0 2 
kX •" co 2 
h 2 - co 2 

He - C0 2 


TABLE 
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0,424 

: 0,803 

" 

1.1*5 

• 0,661 
[ 

61,3 

2,23 

• 

! 0,570 

123,6 

1,23 

| 0,646 

89,1 

2,30 

1 0,750 - 

' 11,7 

2,45 ' 

! 0,720 

MB 
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The eddy mixing coefficient K was chosen within limits 
determined by the theoretical estimates In [10], and then re- 
vised in the process of calculation. 


Several alternate versions were computed: a) with a 

smooth profile for T( /o) and our value of Kj b) with profiles 
of T(f*OandK( ft) from [2]; c) with smooth profiles of T(h) 
and IC(tl) from [2] (auxiliary version to clarify the Influence 
of K on n.(k)). 


Comparison of the profiles (cf. Pig. 1) obtained in ver- 
sions "a" and "b" show that the profiles of FL., (ft.) , obtained /10 

with the smooth profile of r(/0 * and Independent of the altl- 

7 2 ' 

tude of K * 5*10 cm /sec, better describe the experimental data. 

The values of K may decrease sharply at altitudes greater than 
150 -*170 km, i.e., a turbopause may occur, just as on Earth. 

Here, always, ; therefore, the values of K have practic- 

ally no influence on the profile of 

A comparison of "b" and "c" shows the extent to which there 
is little difference between the profiles of h ) when there 
is a large difference in the values of K, especially in the 
case of argon, as follows from the considerations presented 
above. 

5. Let us compare the data obtained from the "Vikings" 
with the data from other experiments and models. According to 
the interpretation of the measurements on the "Vikings" dis- 
cussed above, the most probable temperature of the upper Martian 
thermosphere in the morning (approximately 9 o’clock local time) 
is T = 130 + 30 K, and during the day (approximately 16 o'clock 
local time) is T = 180 + 40 K. If we sum up the experimental 
data from which can be determined (from the altitude scales, 
for hydrogen, carbon dioxide gas and plasma), then, as we showed 
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previously [11, 7], there Is an obvious explicit dependence 
of on the level of solar activity. When solar activity 
ds high, corresponding to a declm.ter solar radio-radiation 
flow for a 10.7 cm wave and P 1Q y * (150 - 200)*10“ 22 w/m 2 Hz, 
the global mean temperature of the upper thermosphere lies 
within the limits of = 350 - 4)0 K; when solar activity 

OO 2$ 2 

is moderately low, corresponding to P^ 0 y = 100 •10" w/rrr 
Ha, 7^ = 270 - 280 K; It is true that this dependency is 
masked at certain periods, obviously, as the result of the 
influence of supplementary thermal heat sources In the thermo- , 
sphere which may be the consequence, on the one hand, of the 
dynamic influence of the lower atmosphere through the upward 
distribution and dissipation into the thermosphere of acoustic- 
gravitational wa v *es, and on the other hand, due to the trans- /II 

mission of energy from the solar wind through the magneto- 
sphere. In every case, the theoretical model [12], in which 
only the basic source of heat due to the solar ultraviolet, the 
infrared heat flow and the elimination of heat by molecular 
conductivity, are taken into account, yields close to the. 
empirical results [11]. Consequently, absorption of solar 
ultraviolet is actually the basic source of heat in the lower 
latitudes of thermosphere in quiet periods. 

At the time of the experiments on "Viking-1" and "Viking- 
2," the level of solar activity was extremely low, correspond- 
ing to F 10< y = 70*10“ 22 w/m 2 Hz [13]. If this is taken into 
consideration, as well as the fact that Mars was close to 
aphelion (which reduces the solar current by approximately 14$ 
in comparison to periods when Mars was at an average distance 
from the sun), then, according to the theoretical model [12], 
we find that in this period, should range from 200 to 220 K. 

We see that nevertheless * is markedly higher than 
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**om the "V’klng" data. How can the still colder thermosphere 
be explained? in our opinion, by the turbulent withdrawal or 
heat downward from the lower thermosphere. Turbulence in the 
lower thermosphere provides, on the one hand, a source of heat 
due to the dissipation of turbulent energy, and on the other, 
it increases the removal of heat due to turbulent thermal con- 
ductivity [7, 14, 15]. To a significant degree, these effects 
are compensated for; however, under certain conditions, one of 
them may predominate. This depends on the magnitude of the 
critical Richardson number which has not yet been sufficiently 
studied, so that contradictory opinions are expressed about 
the predominance of one of the effects mentioned [14, 153. 

As described above, from the "Viking" results in [2], a 
rapid growth in K with altitude is obtained; from (2-4) '10^ 
em 2 /sec at 100 km altitude to (1-4) *10^ cm 2 sec at 170 km 
(Pig. 2). At the same time, according to our interpretation, 
these same data are satisfactorily described in the case of 
I\ at these altitudes — not changing with altitude and such 
that K = 5*10^ cm 2 /sec. 

Let us compare these conclusions winh other data concern- 
ing the eddy mixing coefficient in the upper atmosphere of 
Mars . 


According to our estimate (within the framework of the 
theory of similarity with the help of the Richardson-Obukhov 
law), near the turbopause on Mars, K can change within the limits 
K = (0.9 - 8) *10^ cm 2 /sec [10], Estimates by a similar, but 
somewhat different method in [16], gave values close to our 
lower limit K = 8-10^ cm 2 /sec. Using our estimates [10] and the 
estimates of Golitsin [17], according to which at the surface 
of Mars, K = (0.1 - 4) *10 cmysec, and interpolating between 
them (assuming the linear dependence of fol K on \ ) , we 
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constructed the profile of K( A ^chown in Pic* 2. The homo- 
pause, according to our data, lies at altitudes 120 - 1^0 km. 
Doing these profiles, a model of the composition of the Mar- 
tian atmosphere was constructed [8], the different versions 
of which satisfactorily describe the variations in the con- 
centrations of certain components which were measured experi- 
mentally. On the other hand, the model of Lin and Donahue 
[18], in which the value K = *i-10® cm 2 /sec was used, does not 
depend on altitude and gives reduced values of the concentra- 
tion of 0 and CO near the ionosphere maximum in comparison 
with the experimental values. In [19] the profile of K(^)‘ 
(also represented in Pig. 2) is presented without discussion, 
giving much/ larger values of K than ours. Our attempt to 
calculate a model of the composition [8] makes it possible 

to conclude that these values of K also are overestimated. 

* 

Of course, the profile of K(k ) recommended by us is very 
approximate. Thus, for example, probably at some altitude 
below the homopause (at 150 - 170 km), there exists, just as 
on Earth, a turbopause below which K sharply decreases, since 
on Mars as on Earth, turbulence in the stable lower thermo- 
sphere is maintained due to the energy of the acoustic-gravi- 
tational waves arriving from below and dissipating into the 
lower thermosphere [20-22], However, this has practically no 
influence on the profiles of the concentrations of the com- 
ponents, since here unde' f-ny conditions 

6. Everything discussed above makes it possible to draw 
the following conclusions. 

The convolutions in the altitude profiles of the concen- 
trations of the atmospheric components and the temperature of 
the atmosphere may reflect, to some extent, the influence of 
acoustic-gravitational waves and, to some extent, errors in 


measurements and calculations. In order to isolate the real 
and the artificial waves, a detailed analysis of the methodo- 
logical errors is necessary. The profiles of the temperature 
at the time of the "Viking" experiments may have been smoother 
than the profiles determined in [2] . 

It is possible to determine the eddy mixing coefficient 
from the profiles of the concentrations only when a number of 
conditions are satisfied: stationary distributions of the 

concentrations; a sufficiently large difference between the 
molecular weight of the component used and the mean molecular- 
weight (best of all is to use helium); a sufficiently precise 
determination of the concentrations of all components, includ- 
ing atomic oxygen. The values of the eddy mixing coefficient 
K used in [2] are close to the real values at altitudes of 
100 -'120 km, but at higher altitudes they are too high, and 
the error increases with the altitude. At the time of the 
"Viking" experiments, K = 5*10^ cm^/sec at altitudes of 100 - 
150 km, but at higher altitudes the values of K probably de- 
creased. 

The low temperatures of the thermosphere at the time of 
the "viking" flights (130 + 30 K at 9 o'clock local time, and 
180 + 40 K at 16 o'clock) are explained primarily by the low 
level of solar activity, and to a lesser extent by the fact 
that Mars was close to aphelion, and also possibly to some 
degree by the turbulent removal of heat from the thermosphere. 

Measurements of the vertical profiles for the concentra- 
tions of the atmospheric components, presented with sufficient 
completeness and precision, are some of the most . powerful 
methods for studying 'the structure and the dynamics of a 
planetary atmosphere. 
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CAPTIONS FOR THE FIGURES /l6 

Fig. 1. Height Profiles for the concentrations of components 
of the Martian atmosphere: ; 

\ la) — from "Viking-1"' data; 
lb) — from "Vlking-2" data. 

Experimental points '[1]: 

O - COg, A - Al, H3 - i/ 2 * “ ^2* 

Results of calculations under the following conditions: 

i 7 2 

a) smooth profile of !”(■ K) , . K=5*I0‘cu /c 

b) profiles of T( h ) , K( h ) from [2]; 

c) smooth profiles of from [2 ] ; 


Fig. 2. Altitude profiles of the eddy mixing coefficient K( pV) 
and the coefficient of molecular diffusion fU'N 
in the atmosphere of Mars. 

0' - the values of K from the results of "Viking-1" 
according to [2 ] ; 

A — the values of K from the results of "Viking-2" 
according to [2] ; 

— ■ — — — — the value of K which satisfactorily describes 

the "Viking" data according to th'e present paper; 
■ ■ ■- • — - - — the limits of the possible variations in K 
in the Martian atmosphere according to {8, 10]; 

— . « " ~ — the function K(K ) , recommended in [19]; 
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the function of the coefficients of 
molecular diffusion; 

(a) - "Viking-1", (b) - "Viking-2." 








